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Abstract 

The Casimir energies and pressures for a massless scalar field associated with (5-function poten- 
tials in 1 + 1 and 3 + 1 dimensions are calculated. For parallel plane surfaces, the results are finite, 
coincide with the pressures associated with Dirichlet planes in the limit of strong coupling, and for 
weak coupling do not possess a power-series expansion in 1 + 1 dimension. The relation between 
Casimir energies and Casimir pressures is clarified, and the former are shown to involve surface 
terms. The Casimir energy for a (5-function spherical shell in 3 + 1 dimensions has an expression 
that reduces to the familiar result for a Dirichlet shell in the strong-coupling limit. However, the 
Casimir energy for finite coupling possesses a logarithmic divergence first appearing in third order in 
the weak-coupling expansion, which seems unremovable. The corresponding energies and pressures 
for a derivative of a (5-function potential for the same spherical geometry generalizes the TM con- 
tributions of electrodynamics. Cancellation of divergences can occur between the TE ((5-function) 
and TM (derivative of (5-function) Casimir energies. These results clarify recent discussions in the 
literature. 
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I. INTRODUCTION 



Since the inception of quantum mechanics, divergences associated with zero-point energy 
have caused much confusion. One way to deal with them was to simply define them away. 
This view, however, appears to be untenable, in view of the observable consequence of zero- 
point fluctuations in the Casimir effect, well probed experimentally 0, |^. Calculations 
of such forces, and of the associated energies, are generically plagued with infinities. One 
modern consensus is that Casimir forces between distinct bodies may be unambiguously 
computed, while self-stresses (the very concept of which is only somewhat hazily understood) 
are typically divergent. There are some famous counterexamples: Boyer's result for the 
Casimir energies of a perfectly conducting spherical shell |^, and its generalizations to other 
geometries j^, dimensions i5j, iQ], and fields |7|,J|8|]. Even situations which possess manifestly 
divergent energies, such as a dielectric ball possess unambiguous finite dilute hmits 
1^ 11 1, attributable to van der Waals forces p^ . 

Although these difficulties have been known since at least 1979 [3, ll^ . receritly 
they were rediscovered and reexamined in a series of papers by the MIT group la 17, 3, IS 
HO;!^- Perhaps more heat that light has been generated by some of the recent discussions. 
It is the aim of the present paper to put the discussion on a somewhat clearer footing by 
examining Casimir energies and pressures of massless scalar fields in a (5-function potential 
background. (This is what the MIT group now refer to a the "sharp" limit It is then 

possible to solve the problem exactly, and study how the result depends on the strength of the 
coupling. Although such calculations have been presented by the MIT group [13, E, 2^ 21 
based on the summation of Feynman diagrams, they seem not to have appreciated that 
Casimir energies for such potentials were first computed by the Leipzig group. The first 
calculations with planar 5-function potentials were those of Bordag et al. 22], who found 
equivalent expressions for the Casimir energies given later in Refs. 20, 21j The corresponding 
spherical problem was studied first by Bordag et al. [i^, who found a nonvanishing second 
heat kernel coefficient, indicating that the Casimir energy was divergent in third order 



in the coupling. After a perhaps dubious renormalization, Scandurra [2^ extracted the 



nap 

finite part. Recently, Barton 25|| has carried out related calculations, modeling a FuUerine 
molecule to control and physically interpret the divergences, and examining the TE and TM 
electromagnetic modes, with conclusions not too dissimilar from those of the MIT group. 

Although, therefore, this model seems quite well-studied, it is perhaps worthwhile to re- 
examine it in what I consider the most physically transparent Green's function approach, to 
see if some clarity can be brought to what seems at present a rather confused situation.^ In 
so doing, we shall clarify the discussion of the perturbative expansion, and learn that it is 
only the strong-coupling limit of the spherical Casimir energy that possess a finite self-stress, 
unless cancellations can occur between TE and TM modes (which certainly do occur in the 
strong coupling limit). 

This paper is laid out as follows. In the next section, we find the Casimir pressure 
for a massless scalar interacting with two 5-function potentials in one spatial dimension. 
( Equivalent ly, this is a spherical geometry in one dimension.) The pressure is completely 
finite, but is nonanalytic in the coupling for weak coupling. The Casimir energy receives 



^ Barton refers to my approach as "older methods," but he employs methods of Debye going back to 
early in the previous century, and other classic techniques. I certainly feel in good company if I use the 
propagation functions invented by Green, as well as Debye expansions. 
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contributions from the boundaries (surface terms). The generahzation to (5-function planes 
in three dimensions is immediate, and given in Sec. IIIIl Sec. IIVI presents the corresponding 
calculation for the Casimir energy of a massless scalar interacting with a spherical 5-function 
shell. That resulting expression, in the strong-coupling limit, reduces to the standard one 
for a Dirichlet shell, yielding a finite self-energy j26[. However, for any finite coupling, the 
expression possesses an irremovable logarithmic divergence, which first appears in third-order 



in the weak-coupling expansion pO|, |21|, , although in second order, as noted previously 
[2^, the energy is finite. Section IVl presents the Casimir energy and pressure for a spherical 
derivative of a 5-function potential, which, in the strong coupling limit, corresponds to the 
TM modes of electrodynamics. (The Dirichlet modes computed in Sec. IIVI correspond to 
the TE modes.) Concluding remarks are offered in Sec. IVII 



II. 1 + 1 DIMENSIONS 



We consider a massive scalar field (mass fi) interacting with two 5-function potentials, 
one at a; = and one at a; = a, which has an interaction Lagrange density 



C 



int 



'^-5(a;)0^(x) - ~5{x - a)</)^(x). 



(2.1) 



where we have chosen the coupling constants A and A' to be dimensionless. (But see the 
following.) The Casimir energy for this situation may be computed in terms of the Green's 
function G, 

G{x,x') = i{T(j){x)<p{x')), (2.2) 
which has a time Fourier transform, 



G(x,x')= / ^e-(*-*')^?(x,x». 



which in turn satisfies 



+ + -5{x) + —5{x - a) 
ox^ a a 



g{x, x') = 6{x — x'). 



(2.3) 



(2.4) 



Here = ix^ — uj"^. This equation is easily solved, with the result 



AA^ 

(2Ka) 



-2 cosh k\x — x'l 



A 



for both fields inside, < x,x' < a, while if both field points are outside, a < x, x', 



(2.5a) 



g{x,x 



_}_^-k\x-x'\ _|_ __}_^-K{x+x'-2a) 

2k 2kA 



A 



2Ka 



1 - 



A' 



A' 



2Ka J 2Ka 



1 + 



A 

2Ka 



2Ka 



(2.5b) 
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For X, x' < 0, 



1 



1 



Here, the denominator is 

A 



-k{x+x' —2a) 



\- ' 




- ' ( 


2Ka 




2Ka \ 



2/ta 



2Ka 



X 



2.5c) 
(2.6) 



2/ta J \ 2/ta J (2Ka)^ 
Note that in the strong couphng hmit we recover the famihar results, for example, inside 

sinh Kx< sinh k(x> — a) 



X, X' oo : g{x, x) 



K sinh K,a 



2.7) 



We can now calculate the force on one of the 5-function points by calculating the discon- 
tinuity of the stress tensor, obtained from the Green's function by 



(T^-) = ( d^'d"' - ^g'^'d^d'A ^G{x,x') 



Writing 



we find inside 



2tt ' 



(2.8) 



(2.9) 



1 



txx = -^{uj"^ + d^d-,>)g{x,x'] 
2i 



{2uj^ - /i^ 



4mA 



2/ta 



2Ka 



2Ka 



2Ka 
X' 

2na 



{2Ka) 



X 
2na 



JIkx 



(2.10) 



Let us henceforth simplify the considerations by taking the massless limit, = 0. Then the 
stress tensor just to the left of the point x = a is 



— < 1 + 2 
2i ^ 



2na \ ( 2Ka \ 

+ 1 — + 1 e'"'^ - 1 



A 



A' 



(2.11) 



From this we must subtract the stress just to the right of the point at x = a, obtained from 
Eq. ()2.5b|l . which turns out to be in the massless limit 



x=a+ 



K 

2V 



(2.12) 



which just cancels the 1 in braces in Eq. 1)2.111) . Thus the force on the point x = a due to 
the quantum fluctuations in the scalar field is given by the simple, finite expression 



{Txx) 



x=a+ 47ra^ Jo 
4 



dyy 



1 



{y/X + l){y/X'+l)ey-l 



(2.13) 



-0.10 - 



-0.15 




FIG. 1: Casimir force between two (5- function points having strength A and separated by a distance 
a. 



This reduces to the well-known, Liischer result |27| . I28l | in the limit A, A' 

TT 



OO, 



lim F 

A=A'-^oo 



24a2 



(2.14) 



and for A = A' is plotted in Fig. [T] 

We can also compute the energy density. In this simple massless case, the calculation 
appears identical, because t^x = too (conformal invariance). The energy density is constant 
[Eq. ()2.10p with /i = 0] and subtracting from it the a-independent part that would be 
present if no potential were present, we immediate see that the total energy is = Fa, so 
F = —dE/da. This result differs significantly from that given in Refs. 0,0,^^, which 
is a divergent expression in the massless limit, not transformable into the expression found 
by this naive procedure. However, that result may be easily derived from the following 
expression for the total energy. 



E = J (dr) (r°°) = ^ y (c^r)(a°a'° - V')G(x,x') 



(2,15) 



Integrating over the Green's functions in the three regions, given by Eqs. ()2.5a|) . ()2.5b|) . and 
()2.5c|) . we obtain for A = A', 



E 



dy 



1 



l + y/X 4:na 



1 + 2/{y + A) 

^^^y/A + i)V-r 



(2.16) 



where the first term is regarded as an irrelevant constant (A/a is constant), and the second 
is the same as that given by Eq. (70) of Ref. [3| upon integration by parts. 
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The origin of this discrepancy is the existence of a surface contribution to the energy. 
Because d^T^'^ = 0, we have, for a region V bounded by a surface S, 

= ^ /" (c/r)T°° + / dSiT''\ (2.17) 
dt Jv J s 



Here = d^(j)d^(j), so we conclude that there is an additional contribution to the energy, 



E, = [ dS-VG{x,x'] 

At J 



(2.18a) 
(2.18b) 



where the derivative is taken at the boundaries (here a; = 0, a) in the sense of the outward 
normal from the region in question. When this surface term is taken into account the extra 
terms in Eq. (j2.16p are supplied. The integrated formula (j2.15|) automatically builds in this 
surface contribution, as the implicit surface term in the integration by parts. (These terms 
are slightly unfamiliar because they do not arise in cases of Neumann or Dirichlet boundary 
conditions.) See Fulling ^| for further discussion. 

It is interesting to consider the behavior of the force or energy for small coupling A. It is 
clear that, in fact, Eq. ()2.13|) is not analytic at A = 0. (This reflects an infrared divergence 
in the Feynman diagram calculation.) If we expand out the leading A^ term we are left with 
a divergent integral. A correct asymptotic evaluation leads to the behavior 

\2 1,2 

(ln2A + 7), E -- — (ln2A + 7-l), A ^ 0. (2.19) 



A-Ko^ ' 47ra 

This behavior indeed was anticipated in earlier perturbative analyses. In Ref. the general 
result was given for the Casimir energy for a D dimensional spherical 5-function potential 
(a factor of l/47r was inadvertently omitted) 

o-i-2z>A^ r(i^)r(z^-3/2)r(i-D/2) 

This possesses an infrared divergence as ^ 1: 

E(D=i) ^ ^r(O), (2.21) 
which is consistent with the nonanalytic behavior seen in Eq. ()2.19|) . 



III. PARALLEL PLANES IN 3 + 1 DIMENSIONS 

It is trivial to extract the expression for the Casimir pressure between two 5 function 
planes in three spatial dimensions, where the background lies at a; = and x = a. We 
merely have to insert into the above a transverse momentum transform. 



G{x,x') = / ^e-^"^*-*^ j ^e^'^-(-r)-(7(x,a;';«:), (3.1) 
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where now h? = /i^ + fc^ — cj^. Then g has exactly the same form as in Eqs. ()2.5p . The 
reduced stress tensor is given by, for the massless case, 



T9ix,x') 



(3.2) 



so we immediately see that the attractive pressure on the planes is given by (A = A') 



P 



327r2a4 



dyy^ 



{y/X + iyev - r 



(3.3) 



which coincides with the result given in Refs. j20|, |21| . 

The Casimir energy per unit area again might be expected to be 



S 



1 



dy 



y 



{y/X+iyev 



IP 
3 a ' 



(3.4) 



because then P = —^£- In fact, however, it is straightforward to compute the energy 
density (T°°) is the three regions, z < 0, < z < a, and a < z, and then integrate it over z 
to obtain the energy/area, which differs from Eq. (|3.4p because, now, there exists transverse 
momentum. We also must include the surface term ()2.18aj) . which is of opposite sign, and 
of double magnitude, to the k'^ term. The net extra term is 



1 



487r2a3 



dyy"^ 



1 



{y/X + lyey 



(3.5) 



If we regard A/a as constant (so that the strength of the coupling is independent of the 
separation between the planes) we may drop the first, divergent term here as irrelevant, 
being independent of a, because y = 2Ka, and then the total energy is 



S 



1 



OGvr^a'^ 



dyy^ 



l + 2/{X + y) 
{y/X+iyeV-V 



(3.6) 



which coincides with the massless limit of the energy first found by Bordag et al. |22], and 
given m Refs. 0, |2l|. As noted m Sec. |nj this result may also readily be derived through 
use of ()2.15p . When differentiated with respect to a, Eq. ()3.6|) . with A/a fixed, yields the 
pressure fl3.3|) . 



IV. THREE-DIMENSIONAL SPHERICAL POTENTIAL 

We now carry out the same calculation in three spatial dimensions, with a radially sym- 
metric background 

Ant = -^^5(r-a)02(a;), (4.1) 

which would correspond to a Dirichlet shell in the limit A ^ oo. The time- Fourier trans- 
formed Green's function satisfies the equation (k^ = —u^) 

G(r,r') =5(r-r'). (4.2) 



-V^ + + -5{r - a) 
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We write G in terms of a reduced Green's function 



(4.3) 



Irri 



where gi satisfies 



' ^ ^-^ + K^ + -5{r - a) 



-r" h 

dr dr 



X 
—( 

a 



9i{r,r') = —5{r- r'). 



(4.4) 



We solve this in terms of modified Bessel functions, Iy{x)^ K^{x), where z/ = / + 1/2, which 
satisfy the Wronskian condition 



I'Ax)KJx) - K{x)h[ 



X] 



1 

X 



(4.5) 



We solve Eq. ()4.4|) by requiring continuity of gi at each singularity, r' and a, and the ap- 
propriate discontinuity of the derivative. Inside the sphere we then find (0 < r, r' < a) 



giir, r') 



Krr 



X 



ei[Ka) 



ei[Kry)Si[Kr<) Si[Kr)Si[Kr )— ^ , . , . 

K,a 1 + —si[Ka)ei{Ka) 



(4.6) 



Here we have introduced the modified Riccati-Bessel functions, 



si{x) 



^//+i/2(a;), ei{x) = J'^Ki+y^ix). 




(4.7) 



Note that Eq. ()4.6p reduces to the expected result, vanishing as r — a, in the limit of strong 
coupling: 

1 



lim gi{r, r') 



A^oo 



Krr' 



Kr>)si[Kr^) — -si[Kr)si[Kr ) 

Si[Ka) 



When both points are outside the sphere, r, r' > a, we obtain a similar result: 



gi{r, r') 



Krr 



X 



Si{Ka} 



ei{Kr>}si{Kr<) ei{Kr}ei{Kr ) / n / n 

Ka 1 + —si{Ka)ei{Ka) 



(4.8) 



(4.9) 



which similarly reduces to the expected result as A ^ 00. 

Now we want to get the radial-radial component of the stress tensor to get the pressure 
on the sphere, which is obtained by applying the operator 



drdr' - ^(-9°a'° + V ■ V) ^ Idrdr' - K^ - ^^^^ 



1 

2' 



(4.10) 



to the Green's function, where in the last term we have averaged over the surface of the 
sphere. In this way we find, from the discontinuity of {Trr) across the r = a surface, the net 
stress 

, 00 

^ ^ (4.11) 



27ra2 



1=0 



(e,(x).Ka:))'-^^^^ 
dx 



1 + 



Ae; {x)si{x) 



The same result can be deduced by computing the total energy ()2.15p . The free Green's 
function, the first term in Eqs. ()4.6p or ()4.9|) . evidently makes no significant contribution to 
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the energy, for it gives a term independent of the radius of the sphere, a, so we omit it. The 
remaining radial integrals are simply 



/•oo 1 

y = [(^' + ^(^ + 1)) el + xei4 - x'ef] , (4.12b) 

where all the Bessel functions on the right-hand-sides of these equations are evaluated at x. 
Then using the Wronskian, we find that the Casimir energy is 

^ = -— V(2/ + l)/ dxx — \Yi\l + \lJx)KJx)]. (4.13) 
2'Ka L dx 

i—n U 



If we differentiate with respect to a, with A/a fixed, we immediately recover the force 




This expression, upon integration by parts, coincides with that given by Barton |2^, and 
was first analyzed in detail by Scandurra ^]. It reduces to the well-known expression for 
the Casimir energy of a massless scalar field inside and outside a sphere upon which Dirichlet 
boundary conditions are imposed, that is, that the field must vanish at r = a: 

1 °° r°° d 

lim E = -— y (2/ + 1) / dxx — \n[I^{x)K^{x)], (4.14) 
A^oo 27ia ^-^ In dx 

1=0 ''^ 

because multiplying the argument of the logarithm by a power of x is without effect, corre- 
sponding to a contact term. Details of the evaluation of Eq. ()4.14j) are given in Ref . [i^ . 

The opposite limit is of interest here. The expansion of the logarithm is immediate for 
small A. The first term, of order A, is evidently divergent, but irrelevant, since that may 
be removed by renormalization of the tadpole graph. In contradistinction to the claim of 
Refs. 0, lisi l2fll . I21I 1 , the order A^ term is finite, as claimed in Ref. jl^ . That term is 

E^'"^ = — 5^(2/ + 1) y dxx — [Ii+,/2ix)K^+i/2ix)f. (4.15) 

The sum on / can be carried out using a trick due to Klich js^]: The sum rule 

00 

5^(2/ + l)ei{x)si{y)Pi{cose) = ^e'^ (4.16) 
1=0 ^ 



where p = x'^ + y'^ — 2xy cos 6, is squared, and then integrated over 6, according to 

j dcosePi{cose)Pi,{cose)=6i,^j^. (4.17) 

In this way we learn that 

J2{2l+l)eKx)sKx) = ^ / -e-r (4.18) 

7_n Jo 
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Although this integral is divergent, because we did not integrate by parts in Eq. ()4.15p . that 
divergence does not contribute: 



Ana 



°° , 1 d 
dx -X — 
2 dx 



ix 



dw 



w 



327ra 



(4.19) 



which is exactly the result (4.25) of Ref. [26|, which also follows from Eq. ()2.20j) here. 

However, before we wax too euphoric, we recognize that the order term appears 
logarithmically divergent, just as Refs. |20''] and 2l| claim. This does not signal a breakdown 
in perturbation theory, as the divergence in the D = 1 calculation did. Suppose we subtract 
off the two leading terms. 



1 ^ roo 

E = V(2/ + l) / dx 



X 



d 

dx 



In (1 + XhK,) - \hK, + ^{hK,f 



327ra 



. (4.20) 



To study the behavior of the sum for large values of /, we can use the uniform asymptotic 
expansion (Debye expansion). 



z/ > 1 : h{x)Ky{x) ~ ^ 

2v 



' A{t) Bit) 



(4.21) 



Here x = z/z, and t = 1/v 1 + z"^. The functions A and i?, etc., are polynomials in t. We 
now insert this into Eq. (|4.20j) and expand not in A but in i/; the leading term is 



X' 



00 ^ 



247ra ^ u 
1=0 



dz 



X' 



(1 + ^2)3/2 24™ 



C(i). 



(4.22) 



Although the frequency integral is finite, the angular momentum sum is divergent. The 
appearance here of the divergent ({1) seems to signal an insuperable barrier to extraction 
of a finite Casimir energy for finite A. 

This divergence has been known for many years, and was first calculated explicitly in 
1998 by Bordag et al. [23[, where the second heat kernel coefficient gave 



E 



A3 1 



--, s 



0. 



487ra s ' 

A possible way of dealing with this divergence was advocated in Ref. [j^ ]. 



(4.23) 



V. TM SPHERICAL POTENTIAL 

Of course, the scalar model considered in the previous section is merely a toy model, 
and something analogous to electrodynamics is of far more physical relevance. There are 
good reasons for believing that cancellations occur in general between TE (Dirichlet) and 
TM (Robin) modes. Certainly they do occur in the classic Boyer energy of a perfectly 
conducting spherical shell [3, 31, 32j, and the indications are that such cancellations occur 
even with imperfect boundary conditions jl^. Following the latter reference, let us consider 
the potential 

1 Id 

Cint = 7::Xa-—5{r - a)(j)'^{x). (5.1) 

2 r or 
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In the limit A — > oo this corresponds to TM boundary conditions. The reduced Green's 
function is thus taken to satisfy 



1 d 2 d /(/ + 1) 2 ^ xf 

dr dr r dr 



9i{r,r') = —6{r - r). (5.2) 



At r = r' we have the usual boundary conditions, that gi be continuous, but that its 
derivative be discontinuous, 

_/ , 

n d 



r -r9i 
dr 

while at the surface of the sphere the derivative is continuous, 

r=a+ 



d 
or 



while the function is discontinuous, 



9i 



r=a+ 



0, 



or 



(5.3) 



(5.4a) 



(5.4b) 



It is then easy to find the Green's functions. When both points are inside the sphere. 



r,r' < a : gi{r, r') 



Krr 



Si{Kr<)ei{Kry) 



and when both points are outside the sphere, 

1 



r,r' > a : gi{r, r') 



Krr 



Si{Kr^)ei{Kry) - 



Xnale'ilna)]"^ Si{K,r) si{K,r') 
1 + XKae[{Ka)s'i{Ka) 



A/ta[sJ(/ta)]^e/(/tr)e;(/tr') 
1 + XKae'i{K,a)s'i{K,a) 



(5.5a) 



(5.5b) 



It is easy to see that these supply the appropriate Robin boundary conditions in the A ^ oo 
limit: 

d 



lim —rgi 
A-^o Or 



0. 



(5.6) 



The Casimir energy may be readily obtained from Eq. ()2.15|) . and we find, using the 
integrals ^J^ . 



1 °° f°° d 

E = V(2Z + 1) / dxx—\n[l + \xeUx)sUx)]. 

27ra ^ Jo dx 



(5.7) 



The force may be obtained from this by applying —d/da, and regarding Xa as constant [see 
Eq. ()5.H) ]. or directly, from the Green's function by applying the operator. 



2i 



Vr-Vr' — 



/(/+1) 



9i 



which is the same as that in Eq. (j4.1(jp . except that 

_ 1 



-drr, 



(5.8) 



(5.9) 
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appropriate to TM boundary conditions (see Ref. p], for example). Either way, the total 
stress on the sphere is 

' ^.o, , r,,,. i-iw^w;:, (5.10) 



oo „ 
1=0 



2^2-^ 7o 1 + Xxe[{x)s'i{x) 

The result for the energy (|5.7p is similar, but not identical, to that given by Barton j2^. 
Suppose we now combine the TE and TM Casimir energies, Eqs. (j4.13p and (j5.7p : 

^TE + ^TM ^ L y (2/ + 1) / dxx4- In ffl + A^) (1 + Axe;^;)] . (5.11) 

2na ^-^ Jr. ax i\ x / J 

1=0 ■^'^ 

In the limit A — > oo this reduces to the familiar expression for the perfectly conducting 
spherical shell |3]| : 

li^E = --^yi2l + l) dxx(^ + ^+-^ + -^). (5.12) 

Here we have, as appropriate to the electrodynamic situation, omitted the I = mode. 
This expression yields a finite Casimir energy. What about finite A? In general, it appears 
that there is no chance that the divergence found in the previous section in order A^ can be 
cancelled. But suppose the coupling for the TE and TM modes are different. If A^^A™ = 4, 
a cancellation appears possible. 

Let's illustrate this by retaining only the leading terms in the uniform asymptotic expan- 
sions: (x = uz) 

ei{x)si{x) t i( \ i( \ ^ 

~ — , xe,[x)sAx) ~ , u ^ oo. 5.13 

Then the logarithm appearing in the integral for the energy ()5.1H) is approximately 

The first term here presumably gives no contribution to the energy, because it is independent 
of A upon differentiation, and further we may interpret Yll^o^'^ = Eq. (ICTH) ]. Now if 
we make the above identification of the couplings, 

^ A™ _ 2 

^ - ^ - (5.15) 
all the odd powers of u cancel out, and 

^--^Ep'+i)/ ''-^'■'(1-— )■ (^-is) 

The divergence encountered for the TE mode is thus removed, and the power series is simply 
twice the sum of the even terms there. This will be finite. Presumably, the same is true if 
the subleading terms in the uniform asymptotic expansion are retained. 
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It is interesting to approximately evaluate Eq. ()5.16|) . The integral over z may be easily 
evaluated contour integral, leaving 



oo 

I 

E 



1=0 



This I sum is logarithmic divergent, an artifact of the asymptotic expansion, since we know 
the term is finite. If we expand the square root for small A^/z/^, we see that the (9(A^) 
term vanishes if we interpret the sum as 

oo 

5^z.- = (2^-l)C(.), (5.18) 

1=0 

in terms of the Riemann zeta function. The leading term is (9(A^): 

Y.V^-^- (5,19) 



E ' 

8a ^ v'^ 16a 

1=0 



To recover the correct leading A behavior in ()4.19|l requires the inclusion of the subleading 
y-2n ^grjYis displayed in Eq. 1)4.211) . 

Much faster convergence is achieved if we consider the results with the / = term removed, 
as appropriate for electromagnetic modes. Let's illustrate this for the order A^ TE mode 
(now, for simplicity, write A = A"^^) Then, in place of the energy ()4.19p we have 

^^^ = ^ + — r^sinh2xe-^^ = -f— + —l =-(0.0651061). (5.20) 
327ra 47ra jq ^ \327r 47r / a 

Now the leading term in the uniform asymptotic expansion is no longer zero: 



2na ^ Jo dx\ 8z/2 



= ^E-°(-f) = l^ = -(0-0625), (5.21) 
87ra V 2/ 16a a 

1=1 

which is 4% lower than the exact answer ()5.20j) . The next term in the uniform asymptotic 
expansion is 



a V2048 256 7 a 

which reduces the estimate to 



\2 / o 2 o \ \2 

^ ^ —(0.0027368), (5.22) 



A^ 

^(0) ^ ^(2) ^ _(o.0652368), (5.23) 



13 



which is now 0.2% high. Going out one more term give 

^(4) = [15^(4) _ igi / dzt^ — {7 - U8t^ + 554t^ - 708^^ + 295t*^) 

Sua Jq 16 

59.' +^J!L = _^(0.000158570). (5.24) 



a 524288 a 16328 a 
and the estimate for the energy is now only 0.04% low: 

\2 

^(0) ^ ^(2) ^ ^(4) ^ _(o.06507823). (5.25) 

We could also make similar remarks about the TM contributions. However, evidently 
there are additional subtleties here, so we will defer further discussion for an further publi- 
cation. 



VI. CONCLUSIONS 

In this paper we have repeated some calculations with "sharp" but not necessarily 
"strong" potentials. That is, we have computed Casimir energies in the presence of X6{x — a) 
potentials, in the cases when the delta function lies on two parallel planes (first considered in 
Ref. 01), and when the support of the 6 function is a sphere (first considered in Ref. jisl-li^]). 
We have also considered spherical potentials of the form X6'{r — a)/r. For either spherical 
potential, the approach given here yields finite result in all orders, except the third, provided 
we make the coupling constant identification (|5.15|) in the TM case. That is, the expression 
for the energy possess a logarithmic divergence entirelv associated with the order Feyn- 
man graph. This was rediscovered by Graham et al. 20, 21 1, but obscured by the apparent 



(spurious) divergence they also claimed to find in order A^. The bottom line, however, is 
that these sharp potentials yield a divergent Casimir self-stress. 

The generalizations drawn in Graham et al. papers |2l| are, however, perhaps 



too strong. The fact that the A ^ oo limit of the expression for the energy coincides with 
that for the Dirichlet shell, does not prove that the latter is divergent. It does, however, 
suggest that that idealization does not yield the full result for the energy of a configuration 
defined by a real material boundary. This, of course, is no surprise. It has been recognized 
since at least 1979 0| that constructing a shell from real materials will yield apparent 
divergences as the ideal limit is approached, so for example, a shell of finite thickness made 
of dielectric material will correspond to a divergent Casimir energy. 

So the finite Boyer energy [3J for an ideal sphere results from omitting divergent terms, 
which may or may not have observable consequences. (It may be, of course, that for electro- 
magnetic modes, the divergence found here could cancel, for which we have provided some 
evidence.) However, what is remarkable, and of some significance, is that this finite term 
is unique. For example. Barton has recently exhibited a Buckyball model of a conducting 
spherical shell that possesses various large energy contributions referring to the material 
properties of the shell, but which nevertheless possesses a unique, if subdominant, Boyer 
term of order 1/a |25l |. 

It may be useful to compare this situation with a slightly better understood example, 
the Casimir energy of a dielectric sphere. That is certainly divergent; yet if the divergences 
are isolated in terms that contribute to the volume and surface energies, in order (e — 1)^ a 



14 



unique 1/a coefficient emerges , which may be interpretable as the van der 

Waals energy That coefficient diverges in order (e — 1)^ ji^. This fact seems to bear a 
striking resemblance to the finite Casimir energy found here in order A^, and the divergence 
in the next order. There is also the more than analogous relationship between the finiteness 
of the Casimir energy for a dielectric-diamagnetic ball with e/i = 1, and the finiteness found 
here when ;\TEj^TM _ ^. both cases the divergences separately associated with TE and 
TM modes cancel. 

There are also extremely interesting issues related to surface divergences in the local 
Casimir energy density, which have been discussed recently by Fulling [29|. His ideas likely 
will have bearing on understanding the nature of the divergences encountered in these prob- 
lems. 

Evidently, there is much work to be done in understanding the nature of quantum vacuum 
energy. It would obviously be of great benefit if it would be possible to access these questions 
experimentally. 
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